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Abstract 

We calculate the partition function of a harmonic oscillator with quasi-periodic 
boundary conditions using the zeta-function method. This work generalizes a pre- 
vious one by Gibbons and contains the usual bosonic and fermionic oscillators as 
particular cases. We give an alternative prescription for the analytic extension of the 
generalized Epstein function involved in the calculation of the generalized thermal 
zeta-functions. We also conjecture about the relation of our calculation to anyonic 
systems. 
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Since the generalized zeta-function method for computing determinants was intro- 
duced into theoretical physics calculations PJ, this method has been applied successfully 
in different areas of physics. Examples vary from ordinary quantum and statistical me- 
chanics to quantum field theory where anomalies, fermionic determinants as well as all 
sorts of effective actions can be computed by the generalized zeta-function prescription. 

One of the first and very illustrative examples of this method, published in this Journal 
by Gibbons @], concerns the computation of the partition functions of both the bosonic 
and fermionic oscillators. The application of the zeta-function method for these cases 
demands, the analytical extension of the so called inhomogeneous one-dimensional Epstein 
function (see for instance ref. || for some mathematical details on Epstein functions). 

However, as far as we know, the interpolating system between the bosonic and fermionic 
oscillators has never been treated by the generalized zeta-function method in the liter- 
ature. Recently, this problem has been solved by the authors using the Green function 
method [gj. 

In this letter, we generalize Gibbons work |2j since we compute det(u; 2 — Of) with 
a general boundary condition, which contains the periodic and antiperiodic boundary 
conditions as particular cases, by using the generalized zeta-function prescription. To do 
that, an analytical extension of what is called generalized inhomogeneous one-dimensional 
Epstein function is needed. We solve this problem generalizing the procedure of Ambjorn 
and Wolfran || to the usual Epstein function and obtaining this way a very convenient 
expression for the problem at hand. 

Let us then calculate the partition function for a harmonic oscillator with quasi- 
periodic boundary conditions. For an usual (bosonic) harmonic oscillator and even for 
a fermionic oscillator the corresponding partition functions can be obtained by various 
methods, in particular by a path integral continued to the Euclidean space 

Z(P) = J W x }e~ s , (1) 
where S is the classical action. For the bosonic case we have S = (1/2) J£ x(oo 2 — d^)xdt 



2 



and the integral is taken over periodic functions on the interval (0,/?): 



x(t + /3) = x(t) . (2) 

The fermionic case can be treated on similar grounds substituting the classical (c- 
number) functions x by anticommuting Grassmann variables x and x, which must be 
antiperiodic on the interval (0,/5): 

x(t + (3) = -x(t) ; x(t + P) = -x(t) (3) 

and the path integral measure is [dx] [dx] . 

For the generalized case, however, there is no such a prescription. In order to cir- 
cumvent this difficulty we follow the work of Gibbons who first applied the zeta-function 
method to the calculation of the partition functions for the bosonic and fermionic oscil- 
lators as the following determinants 



Z Bosomc ((3) = det" 1/2 (cu 2 - 8 2 t ) . . (4) 

r p&T%od%c 



and 



rFermionic^ = det +l(^2 _ _ ^ 

CLTht t'p&T' . 

As our starting point, we shall then take as the definition of the partition function for 
the general case, the following determinant 

Z d a ([3) = det°(u; 2 -d?)\ 9 = det°(L) e , (6) 

where a is a parameter related to statistics and 6 to the quasi-periodic boundary condition 
of the chosen oscillator. By quasi-periodic boundary condition we mean that 

x(t + (3) = e w x{t) . (7) 
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The particular cases of the bosonic and fermionic oscillators can be obtained by taking 
a = —1/2, 9 = and a = + 1, 9 = tt, respectively, as we shall show bellow. 

The eigenvalues of the operator L = u 2 — d 2 under the quasi-periodic boundary con- 
dition are: 

Xl = u 2 +(^±^-J ■, (m = 0,±l,±2,...) (8) 
and the corresponding generalized zeta-function is given by 



C\s-L) = £W 



where v = uof3/2ir and we defined 



In 



+oo 

E 



m=—oo 



'9 + Imn 



-2s 



'2tt\~" 9, 
J) 



E 



m=— oo 



(9) 



(10) 



The function ( (s; L) is the generalized thermal zeta-function, which reduces to Gibbons 
cases for particular values of 9. From this function one can calculate the desired partition 
function 



Z e M = exp 



s=0 



Once, the function defined by eq. ( |IT)D is not analytic at the origin we will have to 
make an analytic continuation for T>(s, u, ^) in order to determine the partition function 
Zp(0). In fact, the application of the generalized zeta-function method consists basically 
of the following three main steps: first find the eigenvalues and construct the series ((s) = 
J2 m ^m*; then make an analytical extension for the whole complex s-plane (or at least into 
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s=0 



. For the case at 



a region containing the origin) and finally compute exp j— ^C( s )} 
hand, the function T>(s, u, j-) is in fact related to the one- dimensional generalized Epstein 
function Zf(s; I, a), since 



Zf(s;l,a) = J2[c 2 + (n + 



n=l 



(12) 



for which there are well known analytic continuations M. However, these formulas are 
not easy to handle and we performed a new representation for T>(s, u, j-) which can be 
considered as a generalization of the one given by Ambjorn and Wolfran [||. In our case 
we obtain (see the appendix for details): 



2tt' 



T(s) 



— — — V 4 cos(mi 

m=l 



K s _i(2rmn/) 



(13) 



where K a (x) is the modified Bessel function of the second kind with order a. From 
expression (|l|) it can also be shown that (see the appendix) 



and 



P(0,i/,— ) = 



d ^. 9 s 



In [2 (cosh 2nv — cos ( 



(14) 
(15) 



Substituting (jig ) into flTjJ ), we have: 



exp 



s=0 



exp {a In [2 (cosh 27rz/ — cos 6)]} 
6~ 

2 "~ 2 



4- 



cosh 2 — cos 2 



(16) 



This partition function will reduce naturally to the corresponding one of the usual bosonic 
oscillator if we take a = —1/2 and = 00, 
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2 Bosonic^ = / 2 sinh^J . (17) 

Analogously, for the quadratic fermionic oscillator of Finkelstein and Villasante ]7| we 
have, after choosing a = +1, 9 = ir, that 



^w,X c (/5)=4cosh^. (1 
The linear fermionic oscillator |B| is obtained with the choice L — > L 1 / 2 ; a = +1, 9 = ix 



Z[~ c (P) = 2cosh^. (19) 

In fact, because of the general properties of determinants, the above choice could also 
be a = +1/2, 9 = ir and L unchanged. Curiously, Gibbons was able to find this result for 
the linear fermionic oscillator without a similar rescaling in L. This was possible since he 
considered, from the beginning, the sum in the zeta-function (^) only over positive values 
of m. This naturally falls into our result since this introduces a factor 1/2 for T>(s, u, j-) 
and results in a corresponding power for the determinant, as can be seen, for example, in 
eq. QI3). 

A physical interpretation of our result is also possible. Equation (|16D corresponds to 
the partition function of what we call the anyonic oscillator, in the sense that, besides 
interpolating continuously between the bosonic and fermionic cases it contains these two 
cases as particular choices of the parameters involved. This result coincides with the 
one given by the authors using the Green function method to calculate the ratio of two 
determinants [|J. 

Actually, anyons are expected to live in two space dimensions M, but the generalized 



statistics they are related to, can also be defined in one space dimension [10] as we have 



done here. Furthermore, it has been shown in the literature that anyons confined to the 



lowest Landau level correspond to anyons in one dimension (TT[ . 

As a satisfactory quantization for anyonic systems are not known yet, we believe that 
our result may contribute towards this direction. Besides, we hope that our analytic con- 
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tinuation for T>(s, is, j-) will be of help in discussing other problems, specially in quantum 
field theory as the computation of effective actions at finite temperature, the Casimir 
effect, etc. 

Acknowledgements. The authors were partially supported by CNPq, Brazilian agency. 

Appendix 

Here, we are going to derive the analytic continuation for V(s, z/, given by eq. (|13D, 
in the lines of Ambjorn and Wolfran |J and derive its main properties at s — 0. In fact 
their result corresponds to a particular case of ours, when we take 8 = 0. 

Starting from eq. (|T0|) and using the integral representation for the Gamma function 





T(s)A- 

which is valid for Re r > we find 



dr r s - l e~ AT 



(20) 



2tt' 



E 



m=— oo 



2tt 



+ m 



I POO 

— dx X s ' 1 exp 

1 f°° r + °° f / 

— / dx x s ~ l exp I —v 2 x \ exp < — 



2tt 



X 



+ m \ x 



Now, using the Poisson summation formula [|13 | 



+ 0O 

E e 

n=—oo 



-(z+n) 2 ^ 

and identifying m = n, z = 0/2tt, x = tt/t we can write 



r E e " 

n=— oo 



n 2 7TT—2inzn 



(21) 



(22) 
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9 1 f°° r -i pK + °° { 

2tt r(s) Jo - 1 V x m ±^ OD [ 



,7T 



m im^ 



+oo 

5] e" imy / f/.r.r^-- J ( 



, OO . , , 2 

im# / j„ „«— 4— i/ as— m — 



(23) 



The integral in the last line of the above equation can be written in terms of a modified 
Bessel function of the second kind, K a (y) 112] 



(24) 



which is valid for Re/3 > and Rcy > 0. Then, splitting the sum over m from — oo 
to +oo into the sums over negative values, m = and positive values and identifying 
a = s — 1/2, (3 = m 2 iT 2 and 7 = u 2 , we have 



^ 



(25) 



where we defined 



G(s) = v^<|4^ + 4 £ cos ( m6 

m=l 



K s _i(2m7ru) 



(26) 



± ' . ' \ J/ 

which is the result expressed in eq. fljjl) . 

Now, we are going to compute P(0, z/, j-) as well as &D(0, v, j-)/ds. Noting that G(s 
is an entirely analytical function and that T(s) diverges as s — > 0, we see that 



P(0, 1/, = . 

Next, we turn to the derivative of the zeta-function ( e (s; L) at the origin, 



(27) 



9s 



C*(s;£) 



s=0 



d 9 , 

— Pis, za — 
9s V ' '2tt' , 



s=0 



(28) 
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where we already used the fact that P(0, v, j-) = 0. Using eq.(|25|) we have 



d ^, e , 



r*(s 



G(s) 



+ 



r s) 



-G'(s) 



s=0 



and the second term of the r.h.s. of the above expression vanishes too. Noting that 



T'(s) 



lim ^ , 

s->o T 2 (s) 

and using the fact that r(— |) = —2y/n, we have 



d 9 , 

— Pis, v, — 
ds v ' '2tt' 



G(0) 



/ J/ \ 2 

-27Ti/ + 4 ^ cos(m0) I — J (2m7rz/) 
m=i Vm/ 



Then, using the well known expression for K_i (z), namely, [12 



and the result 



we finally obtain 



KAz) 



OO 1 

n 

n=l n 



7tv — ln[2 sinh(7r^)] , 



9 <n( 6 



—2-kv + 2 cos(m6) — e 



s=0 



m=l 
-oo ^ 



-2nu + ~ e 



in 



-(2irv—i9)m 



m=l 



m 



-2m-KU 



+ E 

m=l 



-(27w+j0)m 



rn 



In 



2 sinh 7tz/ + 



i9" 



= — In [2 (cosh(27rz/) — cos i 
which is precisely the result shown in eq. (|T5D . 
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